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1 Introduction 



The Calogero-Moser-Sutherland-type many-particle models are intensively studied integrable systems 
with deep connections to many important branches of mathematics and physics. As a classical Hamil- 
tonian system, the hyperbolic BC n Sutherland model is defined on the cotangent bundle of the open 
subset 

C = {q= ( Ql , . . . ,q n ) £ R n | Ql > . . . > q n > 0} C R n . (1.1) 
For convenience we identify the phase space T*c with the product manifold 

V s = C xR n = {(q,p)\q € c, pel"}, (1.2) 

endowed with the standard symplectic form 



^2<iq c Adp c . (1.3) 



c=l 

The dynamics is governed by the interacting many-body Hamiltonian 



HS = 2^ + ^ (^mhfe) + sinh 2 (2, c ) J 

9 2 



E 



c=l 

~2 J2 



+ 



i<a<b<n vSinh 2 (g a -%) smh 2 (q a + q b ) 



(1.4) 



with coupling constants g 2 ,gf,g 2 £ M satisfying g 2 > 0, g\ > and g 2 > —\g\- 

By applying the projection method on the geodesic system of the non-compact Riemannian sym- 
metric space SU(n + l,n)/S(U(n + 1) x U(n)), Olshanetsky and Perelomov constructed a Lax rep- 
resentation of the BC n Sutherland dynamics and analyzed the issue of solvability as well, but only 
under the restrictive assumption g\ — 2g 2 + \/2gg2 = (for details see e.g. pQ, [2], [3]). As the algebraic 
methods prevailed, the Lax representation of the dynamics was soon established for arbitrary values 
of the coupling constants (see e.g. [1], [S], [B], [7]). Somewhat surprisingly, the symplectic reduction 
derivation of the BC n Sutherland model with three independent coupling constants is only a relatively 
recent development [5]. Besides providing a nice geometric picture and an efficient solution algorithm, 
the symplectic reduction approach has also allowed us to construct action-angle variables to the BC n 
Sutherland model and to establish its duality with the rational BC n Ruijsenaars-Schneider-van Diejen 
system (see [9]). Sticking to the powerful machinery of symplectic reduction, in this paper we con- 
struct a dynamical r-matrix for the BC n Sutherland model. By accomplishing this task we generalize 
the results of Avan, Babelon and Talon on the r-matrix structure of the hyperbolic C n Sutherland 
model, which appeared in their paper [TO] , 

The rest of the paper is organized as follows. In order to keep the presentation self-contained, in the 
next section we provide a brief account on the group theoretic and symplectic geometric background 
underlying the symplectic reduction derivation of the hyperbolic BC n Sutherland model. Built upon 
the reduction approach outlined in Section [21 in Section [3] we construct a ^-dependent dynamical 
r-matrix for the most general hyperbolic BC n Sutherland model with three independent coupling 
constants. The new results are summarized concisely in Theorems 2] and [5j Subsequently, in Section 
[31 we offer a short discussion on possible applications and related open problems. Finally, some 
auxiliary material on the Lie algebra u(n, n) can be found in an appendix. 
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2 Preliminaries 

In this section we review the symplectic reduction derivation of the hyperbolic BC n Sutherland model. 
For convenience, we closely follow the ideas and conventions presented in [9]. 

2.1 Group theoretic background 

Take an arbitrary positive integer n S N, let N = 2n, and consider the N x N matrix 



On In 
In On 



(2.1) 



The matrix Lie group 

G = {yeGL(N,C)\y*Cy = C} (2.2) 
provides an appropriate model of the real reductive Lie group U(n,n). Its Lie algebra 

= u (u, n) = {Y 6 Q l(N, C) | Y*C + CY = 0} (2.3) 
comes naturally equipped with the Ad-invariant symmetric bilinear form 

<,): fl Xfl->M, (Y,Y)^(Y,Y) = tr(YY). (2.4) 

The fixed-point set of the Cartan involution @(y) = (y -1 )* {y E G) can be identified as 

K = {y£G\e(y) = y}^U(n)xU(n), (2.5) 

meanwhile the corresponding Lie algebra involution 9(Y) = —Y* (Y E g) naturally induces the Cartan 
decomposition g = t © p with the eigenspaces 

e = ker(0-ld g ) and p = ker{9 + Id g ). (2.6) 

That is, each Y € q can be uniquely decomposed as Y = Y + + Y_ with Y + € t and Y_ G p. Note that 
the bilinear form (|2.4|) is negative definite on the subalgebra t, whereas it is positive definite on the 
complementary subspace p. 

Now, with each real n-tuple q = (qi, . . . , q n ) E W 1 we associate the N x N diagonal matrix 

Q = diag(gi,. . .,q n ,-qi, . . .,-q n ) G P- (2.7) 

Clearly the subset o = {Q £ p\q £ M. n } is a maximal Abelian subspace in p, which can be naturally 
identified with R n . Under the adjoint action, the centralizer of a in K is the subgroup 

M = Z K (a) = {diag(e iXl ,...,e ix ",e iXl ,...,e ix ")|xi,...,Xn CK (2.8) 
with Lie algebra 

m = {diag(ixi, • • • , iXn, i-Xi, ■ ■ ■ , ' l Xn) \ Xi, • • • , Xn G M} C t (2.9) 

Let (respectively m^) denote the set of the off-diagonal elements of p (respectively t); then with 
respect to the bilinear form (|2.4p we have the refined orthogonal decomposition 

= mem- 1 © a© cr 1 . (2.10) 
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Practically, each Lie algebra element Y £ q can be decomposed as 

y = Y m + y m x + Y a + y a x (2.11) 

with unique components belonging to the subspaces indicated by the subscripts. 

Notice that for each q G W 1 the operator ad<g = [Q, •] G fll(fl) leaves the subspace m J -©a J - invariant; 
therefore the restricted operator 

adg = ad Q | m x ea x G g^m 1 ® o x ) (2.12) 

is well defined. Recall that the regular part of a consists of those diagonal matrices Q, for which the 
linear operator ad<g is invertible. Clearly the standard Weyl chamber 

{Q G a|gi > ... > q n > 0} (2.13) 

is a connected component of the regular part of a. For simplicity, in the rest of the paper we shall 
identify this Weyl chamber with the configuration space c (jl.ip . 

In deriving the Sutherland model from symplectic reduction, the so-called KAK decomposition of 
G plays a crucial role. It basically says that the map 

R n x K x K B (q, k L , k R ) t-> k L e Q k^ G G (2.14) 

is onto. Let G reg denote the image of c x K x K under the above map. As is known, the subset of 
regular elements, Greg, is an open and dense submanifold of G. Moreover, the smooth map 

7r: c x K x K -» G rcg , (q,kL,kn) i-> k L e Q k^ (2.15) 

is a principal M-bundle in a natural manner. Consequently, we arrive at the natural identification 
G reg — c x [K x K)/M*, where M* stands for the diagonal embedding of M (|2.8j) into the product 
Lie group K x K. That is, M* consists of the pairs (m, m) G K x K with m G M. 

We conclude this subsection with a brief excursion on certain adjoint orbit of fi, which is at the 
heart of the symplectic reduction derivation of the BC n Sutherland model. For this, let us consider 
the following set of column vectors 

S = {V €C N \CV + V = 0,V*V = N}, (2.16) 

which can be seen as a sphere of real dimension 2n — 1. With each vector V G S we associate the 
matrix 

£{V) = ifJ,(VV* - 1 N ) + - v)C G t, (2.17) 

where |i,i/£l \ {0} are arbitrary non-zero real parameters. Let us now introduce the distinguished 
column vector E G S with components 

E a = -E n+a = l (o€N n = {l,...,n}). (2.18) 

Notice that the adjoint orbit in t passing through the element £(E) has the form 

= 0(£(E)) = {£(V)et\V £S}. (2.19) 

More precisely, with the free action U(l) x S 3 (e 1 ^, V) i— )■ e l ^V G S, the map 

C:S^O, V^£{V) (2.20) 

is a principal [/(l)-bundle. Therefore the identification O = S/U(l) is immediate. 
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2.2 The Sutherland model from symplectic reduction 

We continue with a short review on the symplectic structure of the cotangent bundle of G. For conve- 
nience, we trivialize this bundle by left translations. Therefore, making use of the linear isomorphism 
g* = g induced by (|2.4p . we can think of the product manifold V = G x g as an appropriate model of 
T*G. At each point (y, Y) G V the canonical symplectic form uj G £l 2 (P) can be written as 



(y,Y) 



(Ay © AY, Sy © 5Y) = (y^Ay, 5Y) - (y^Sy, AY) + ([y^Ay, y^Sy^Y), (2.21) 



where Ay® AY and 5y®5Y are arbitrary elements belonging to the tangent space T y GOg = T^ yY )P. 
Turning to the adjoint orbit O (|2.19|) . remember that it also carries a natural symplectic structure 
induced by the Kirillov-Kostant-Souriau symplectic form uj° G Vt 2 (0). Let us keep in mind that at 
each point p G O it takes the form 

u>%([X,p],[Z,p]) = (p,[X,Z}), (2.22) 

where [X, p], [Z, p] £ T p O are arbitrary tangent vectors with X, Z £ t. Now, motivated by the standard 
'shifting trick' of symplectic reduction, we introduce the product symplectic manifold 

(V ext ,u ext ) = (VxO,uj + uj ). (2.23) 

For an arbitrary function F G C°°('P ext ), at each u = (y,Y, p) G p ext ) we define its gradients 

V G F(u) G g, V s F(u) G g, V°F(u) G T p O C 6, (2.24) 

by requiring 

(dF) u (5y ®5Ye [X, p\) = {V G F(u), y'Hy) + {V s F(u), SY) + (V°F(u),X) (2.25) 

for all 5y G T V G, 5Y G g and X G t. By combining the definition Xp joj cxt = dF with the above 
formula, for the Hamiltonian vector field Xp G £("P cxt ) we find 

(X F ) U = (yV*F(u)) y © ([Y, V s F(u)} - V G F(u)) Y © (-V°F(u)) p . (2.26) 

Consequently, from the definition {F, H} ext = u cxt (Xp, Xh) it is immediate that the Poisson bracket 
of any pair of functions F,H£ C°° (V ext ) takes the form 

{F,Hr x \u) = (V G F( U ),V»H(u)) - {V G H(u),V*F(u)) 

- {[V s F(u),V s H(u)],Y) + u p 9 (V°F(u),V°H(u)). [ ' ' 

Inspired by the KAK decomposition of G ()2.14j) . let us observe that the map 

<D cxt : (K x K) x V cxt ^V cxt , ((k L ,k R ),(y,Y,p)) ^ (k L yk R \k R Yk~ R \k L pkl l ) (2.28) 

is a symplectic left action of K x K on J> ext } admitting a K x i^-equivariant momentum map 

jext . pext ^ t ^ ^ y- ^ ^ ((yYy- 1 ) + + p) © (-Y+ - KiC) (2.29) 

for all s6l. As is known (see [8], [9]), the phase space of the hyperbolic BC n Sutherland model can 
be derived by reducing the extended phase space J> ext at the zero value of the momentum map J ext . 
In the following we briefly summarize the main steps of the reduction. 
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First, let £o denote the set of those points u of the extended phase space (|2.23p . for which we have 
J ext (u) = 0. Note that the level set £o turns out to be an embedded submanifold of ~p ext . To analyze 
its finer structure, we introduce the Lax matrix 

L : V s -> , L(q, p) = L p (q, p) - K iC, (2.30) 

where L v (q,p) £ p is an Hermitian matrix having the block matrix structure 



A B ' 

-B -A 



(2.31) 



More precisely, the entries of the n x n matrices A and B are defined by the formulae 



A ZIHl a -„ n [ _t n - " + KCOsh ( 2 Qc) , , 

smh(g a - g 6 ) smh(g a + g 6 ) smh(2g c ) 

where a,b,c £ N n , a / b. We also need the manifold M s = V s x {K x K)/U (1)*, where f/(l)* denotes 
the diagonal embedding of C/(l) into K x K. Now, one can verify that the map 

T 5 : .M 5 ^P ex \ (q,P,(VL,VR)U(l)*) ^ (VLe Q VR 1 ,mL(q,p)VR 1 ,VLaE)v L 1 ) (2-33) 

is an injective immersion with image Y S (A4 S ) = £q. Since the manifolds A4 S and £o are of the same 
dimension, this observation leads to the identification £o = Ai s . 

Second, by examining the (residual) action of K x K on the model space M of the level set £o, 
it is immediate that the base manifold of the trivial principal (K x K)/U{\) ^-bundle 

tt S :M S ^V S , (q,P,(r}L,VR)U(l)*)^ (q,p) (2.34) 

provides a convenient model for the reduced symplectic manifold. That is, we end up with the natural 
identifications 

V ext // (K x JQ = M S /(K xK)^ V s . (2.35) 

Making use of the defining relationship (ir s )*u} vcd = (T s )*oj cxt , for the reduced symplectic form we 
find immediately that w red = 2lo s with the canonical symplectic form u (jl.3p . Consequently, for the 
reduced Poisson bracket we obtain 

{.,.} s = 2{.,-Y ed . (2.36) 
Finally, let us consider the K x if-invariant quadratic Hamiltonian 

F 2 (y,Y,p) = \(Y,Y) = \tr(Y 2 ) ((y,Y,p) 6 P ext ). (2.37) 

It is clear that F% generates the 'free' geodesic motion on the group manifold G. Due to its invariance, 
it survives the reduction and the corresponding reduced Hamiltonian coincides with the Hamiltonian 
of the Sutherland model (jl.4[) with coupling constants 

g 2 = p 2 , g\ = \vK, g\= X -{v-K) 2 . (2.38) 

Just now can we really appreciate the inclusion of the innocent looking K-dependent central element 
KiC into the momentum map J cxt (|2T29L Indeed, by special izing the parameters (/i, z^, k) appropri- 
ately, from the proposed reduction picture we can recover the most general hyperbolic BC n Sutherland 
model with three independent coupling constants. 
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3 Construction of the r-matrix 



Since the eigenvalues of L (|2.3U|) are in involution (see [3]), we know from general principles that 
the Lax matrix obeys an r-matrix Poisson bracket. As is known (see e.g. [H]), there is a general, 
purely algebraic approach to find an explicit formula for the r-matrix. However, we rather follow the 
symplectic reduction approach put forward by Avan, Babelon and Talon in [10]. It is worth mentioning 
that this geometric approach was later generalized and systematically exploited in [12], leading to a 
uniform treatment of the classical r-matrix structure for various integrable systems. 

3.1 Local extensions of the Lax matrix 

Take an arbitrary point (q,p) of V s and keep it fixed. Notice that the point 

z = (q,p,(l N ,l N )U(l)*)£M S (3.1) 

projects onto (q,p), i.e. tt s (zq) = (q,p). Now, pick an arbitrary function / G C°°(V S ). We say that a 
smooth function 

f-.U—t-R (3.2) 
defined on some open neighborhood U C *p ext of point 

u = T s (z ) = {e Q ,L{q,p)^{E)) G P ext (3.3) 
is a local extension of f around uq, if 

/oT 5 | (TS) _ 1(c/) = /o^| (TS) _ 1([/) . (3.4) 

As is known, this special class of local extensions can be used effectively to compute reduced Poisson 
brackets by evaluating certain 'unreduced' Poisson brackets. More precisely, if f,g G C°°(U) are 
arbitrary local extensions of functions /, g G C°°(V S ) around uq in the sense of (|3.4p . then 

{f,gy cd (Q,p) = {f,gr t (u ). (3.5) 

In particular, in the following we shall make use of the above formula[ll to find an explicit expression 
for the r-matrix of the Sutherland model. The auxiliary functions defined below play an important 
role in constructing local extensions of the Lax matrix L (|2.30p . 

We start with the study of the smooth principal M-bundle ir (|2.15|) induced by the KAK decom- 
position of G. Since 7r(g, ljv, ljv) = e , there is a smooth local section 

G 3 y ^ (a c (y),a L (y),a R (y)) ecxKxK (3.6) 

of 7T, defined on some small open neighborhood G C G7 rog of e**, such that 

{a c {e Q ),a L {e%cx R {e Q )) = {q,l N ,l N ). (3.7) 

Besides the above normalization, we may impose certain conditions on the derivative of section (|3.6p 
at , too. Notice that the tangent space of c x K x K at (q, In, In) can be identified as 

T (g ,i N ,i N) (c x K x K) = T q c ® T lN K ® T 1n K * R n ® * t, (3.8) 

Consistently with the Dirac bracket, generally there are also some correction terms on the right hand side of (13.51 1. 
However, since we reduce at the zero value of the (equivariant) momentum map, and since by (|3.4fl our local extensions 
are (locally) K x Jf-invariant on the level set £o, these corrections cancel. For details see e.g. Chapter 14 in |13j . 
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in which the vertical subspace of it takes the form 

ker((dvr) (giljv>ljv) ) = {Oel©l£l n et©t|Iem}^m. (3.9) 

Since K n © © {! is clearly a complementary subspace of the vertical subspace, we may assume that 
at point the derivative of (j3.6j) maps into the complementary 'horizontal' subspace, i.e. 

ran((d(a c , a L , <r R )) eQ ) = R n © m ± © t (3.10) 

That is, we may assume that 

ran((do-L) e g) = m" 1 . (3-11) 
In the following we will need an explicit formula for the derivative of (|3.6p at point . 
Lemma 1. Under assumption H3.11\) . for each tangent vector 5Y G g = T\ N G we have 

(da L ) eQ (e Q 5Y) = - sinh^d Q )-\SY) a ±, (3.12) 

(da R ) eQ (e Q 5Y) = -(5Y) + - coth(ad Q )(5F) a x . (3.13) 

Proof. For simplicity, let us introduce the shorthand notations 

5 c = (da c ) eQ (e Q 5Y)GR n , 5 L = {da L ) eQ {e Q 5Y) G m ± , 5 R = {da R ) eQ {e Q 5Y) G 6, (3.14) 
and define D c = diag(5 c , — 8 C ) G o. Since G is open, for small values of \t\ we have e®e tSY G G, whence 

e®e tSY = (Ti (e Q e t ' 5y )e dia ^( eQe ^)'-^( eQe " y ))^(e Q e t5y )- 1 . (3.15) 
By taking the derivative of the above equation at t = 0, we obtain 

5Y = cosh(ad<g)5L - sinh(ad<2)c)L + D c - 5 R . (3.16) 
It follows that (5Y) a = D c and (5Y) m = — (5 R ) m , meanwhile for the off-diagonal components we get 

(SY) a ± = -sinh(adQ)5L, (SY) m ± = cosh(adg)5 L - (S R ) m ±. (3-17) 

By solving this linear system for L> c , 5l and 5 R , the lemma follows. □ 
To proceed further, let us note that 

8 = {VeS\V 1 ^0,...,V n ^0] (3.18) 
is an open and dense submanifold of S (j2. 16[) . which contains E (|2.18[) . The map 

r:S^M, ^diag^,...,^,^,...,^) (3,9) 

defined on S is smooth, satisfying t(E) = In- Utilizing the natural identification 

T E S ^ {5V G | C5V + 5V = 0, (5V)*E + E*5V = 0}, (3.20) 
for the action of the derivative of r on each tangent vector 5V G T E S we find 

(dT) E (5V) = idiag (Im(OnOi), . . . , Im(((jy) n ), Im((<N0i)> . . . ,lm((SV) n )) G m. (3.21) 



Turning to the principal [/(l)-bundle £ ()2.20p . notice that E projects onto £(E). Therefore, we 
can find a smooth local section 

6 3 p ■-> V(p) e 5 (3.22) 

of £, defined on some open neighborhood O C of such that V(£(i£)) = E. Moreover, by 

'shrinking' O if necessary, we may assume that 

V(p)£S (VpG0). (3.23) 

In order to fix the range of the derivative of section V at point £(E), notice that the map 

T E S x T E S B (8V, 5W) i — y Re ({5V)*6W) G R (3.24) 

is an inner product on the tangent space T E S (|3.20p . Therefore, we may assume that the derivative 
operator (dVW^) maps T^ e ^O into the orthogonal complement of the vertical subspace 

ker((d£)i?) = Mi£ = {xiE G T E S\x G R}. (3.25) 

This requirement amounts to the constraint 

(SV)*E = E*5V (ySV G ran((dV)£(£))). (3.26) 

Remembering the local section (|3.6p . let us consider the smooth map 

7:GxO->0, (y, p) ^ ( J L (y)- 1 p^L(y)- (3-27) 

Notice that j(e^,^(E)) = £(E) G O. Therefore, there are some open neighborhoods G C G of e^, 
and C 6 of f (£), such that 

T(W)£0 (V(w)eGxO). (3.28) 
Having equipped with r, V and 7, now we define their composition 

m:Gx6^M, (y,p)^ m{y, p) = r(V( 7 (y, p))), (3.29) 

which is a smooth map satisfying m(e < *,£(i£)) = Ijv- To compute the derivatives of the matrix entries 
of the diagonal matrix 

m(y,p) = diag(mi(y,p), . . . ,m n (y,p),m 1 (y,p), . . . ,m n (y,p)) (3.30) 

at point (e® ,£(E)), for each c G N n we introduce the column vector F c G with components 

(F c ) a = -(F c ) n+a = 6 c>a (a G N n ), (3.31) 

together with the Lie algebra element 



i ( F C E* + EF* - ^EE* ) G f . ( 3.32) 



Lemma 2. Ta&e arbitrary vectors SY G g and Z G t; i/ien for each c G N n we /ioue 

(dm c ) {eQ ^ E)) (e Q 5Y(B[Z,C(E)]) = (s c , Z + sinh(S Q )- 1 (^) x) • (3.33) 



Proof. From (|3.29p it is clear that 

(dm) (eQ>m) ^6Y 8 [Z,£(E)]) = (dr) E (dV) aE) (d 1 ) {eQ ^ E)) (e^5Ye[Z,aE)})- (3-34) 
Upon introducing the shorthand notation 

X = Z + sinh(ad Q )" 1 (5y) a ± G C, (3.35) 
from ()3.27p and Lemma [1] it is immediate that 

m (eQ ^ E)) {e Q 5Y ®[Z^(E)}) = G T^ E) 0. (3.36) 

Next, let us consider the tangent vector 

5V = (dV) e(£;) (d7)( e Q,^ ) )(e Q <5y © [Z, £(£)]) G T £ S. (3.37) 
Since V ()3.22p is a local section, we have £ o V = Id^, which entails 

(dO B *V = (d£) £ (dV)^)[X,£(£)] = (3.38) 

On the other hand, notice that the vector XE belongs to the tangent space T E S (|3.20p . Moreover, 
we find easily that 

{dOE{XE) = [X,aE)]. (3.39) 

From the last two equations we conclude that W — XE belongs to the vertical subspace (|3.25p . whence 
5V = XE + x\E with some x G M. However, the value of x is uniquely determined by (|3.26p . from 
where we infer that 

Tp* V Tp 

5V = XE- ^^E. (3.40) 
Now, by combining (j3.34j) and (|3.37p . we see that 

(dm) (eQ ^ E)) (e Q 5Y(B[Z,aE)]) = {dr) E 8V. (3.41) 
Therefore, recalling (|3.2ip and (|3.3U|) . for each c G N n we can write 

(dm c \ eQ>m (e^5Y © [Z^{E)\) = iM(«0c) = ilm((^e) - ^ X ^\ (3.42) 

Utilizing (|3.31j) . notice that 

iIm((XE) c ) = hlm(F*XE) = ^{F*XE + E*XF C ) = jM(i(F c E* + EF*)X). (3.43) 

Plugging this formula into (|3.42p . the lemma follows. □ 

Now we are in a position to construct an appropriate local extension of the Lax operator L (|2.30p . 
For this, let us consider the open subset 

p ext = G x g x 6 C V ext , (3.44) 

which clearly contains the reference point uq (|3.3p . Recalling ()3.6p and ()3.29p . on p 6 ** we define the 
smooth function 

tp-.-pBt^K, (y,Y,p)^ip(y,Y,p) = a R (y)m(y,p). (3.45) 
Notice that <p(uo) = In- Finally, let us consider the locally defined smooth function 

L:V ext ^Q, (y,Y,p)^L(y,Y,p)=^(y,Y,p)- 1 Y^(y,Y,p). (3.46) 
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Lemma 3. The Q-valued function L \3.J^6^ , defined on a small open neighborhood of the reference 
point iio h 6 3. 6 J\) , is a local extension of L /12.3U\) . More precisely, we have 

° '^ 5 '!(T s )- 1 (7' cxt ) = L ° \tys^-i(pexty (3.47) 

Proof. Take an arbitrary point z = (q,p, fj R )U(l)*) G (Y 5 ')" 1 ^ 6 **), and let 

(y,Y,p) = T S (S) = (fj L eQfj-\fj R L(q,p)fj-\fj L Z(E)fjl l ) € V cxt . (3.48) 

By applying the local section a (|3.6|) on the Lie group element y £ G C G, we see that 

V = cT L {y)e^^-^a R {y)-\ (3.49) 

Recalling ir ()2.15p . it is immediate that q = cr c (y) and (fjL, fjjijM,, = (o-L(y),cr R (y))M*. Therefore, we 
have 

fjL = o- L (y)V and fj R = o- R (y)V (3.50) 

with some diagonal matrix T> £ M. 

Next, we inspect the Lie group element m(y,p) = t(V(7(?/, p))) £ M. Remembering (|3.27|) . notice 
that 

7 (y, p) = a L {y)- x pa L {y) = vUvV^Ltfflrfc l a L {y) = V^E^ 1 = £(VE). (3.51) 
Since V (|3.22|) is a local section of £ (|2.20p . we find that 

£(V( 7 (y, p))) = {i o V){i{VE)) = i(VE). (3.52) 

Therefore, there is some i> e R, such that V(>y(y,p)) = e^VE. Recalling ([3TT9D . it follows that 

m(y, p) = T {e^VE) = e^V. (3.53) 

Due to relationships (|3.5Up and (|3.53|) we observe that 

<f(y, Y, p) = a R (y)m(y, p) = e l ^f] R . (3.54) 

Therefore, recalling (|3.46p and ()3.48p . we find immediately that 

L(T s (z)) = <p(y, Y, py'Y^y, Y, p) = L(q,p) = L(tt s (2)). (3.55) 

Since z is an arbitrary element of (T' s ) _1 ('P ext ), the lemma follows. □ 

3.2 Computing the r-matrix 

Let us choose some dual bases {7a}, {T a } in g, i.e. (T A ,T R ) = S A , and consider the function 

V ext 3 (y, Y, p) ' y Y(y, Y, p) = Y G , (3.56) 

together with its components 

V ext 3 (y,Y,p) ' y Y A (y,Y,p) = (T A ,Y) G M, (3.57) 
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defined on the extended phase space. Notice that Y = J2a yAt a- 

Since Y depends only on variable 

Y, its only nontrivial gradient (|2.24|) is 

V S Y A = T A , (3.58) 

whence from ([237]) we obtain that {Y A . Y B } ext = —([T A ,T B ],Y). As usual in the theory of integrable 
systems, these relationships can be succinctly rewritten in the standard St Petersburg tensoriaH 
notation. Indeed, upon introducing the quadratic Casimir 

Oi2 = ^21 = ^T A ®T A £ 5 ®g, (3.59) 

A 

we find easily that 

{Y 1 ,Y 2 } cxt = ^{1^, Y B } cxt T A ®T B = [-0 12 /2, Fi] - [-0 2 i/2, Y 2 ]. (3.60) 

A,B 

Now, from (j3.46p and ()3.56j) we see that L can be obtained from Y by the gauge transformation 

L = ip- 1 Yip. (3.61) 
It readily follows^] that {L\, L 2 } ext = [r"i 2 , L\] — [f 2 i, L 2 ] with the transformed r-matrix 



fl2 = f 1 



-inia - {^ 1) r 2 } ext vr 1 + 5 [{^ 2 } cx VrV 2 ~\^ 2 ]) (pn^. (3.62) 



Recalling (|2.36p and (|3.5|) . it is thus immediate that for the Lax matrix L (|2.3Up we have 

{L ll L 2 } S (q,p) = [r 12 (q,p),Li(q,p)} - [r 21 {q,p), L 2 (q,p)] (3.63) 

with the r-matrix 

r 12 (q,p) = 2r~ 12 (n ) = -fii 2 - 2{<p u Y 2 } cxt (u ) + [fa, </> 2 } ext (u ), L 2 {q,p)] . (3.64) 

In the following we use extensively the special basis of g introduced in the appendix. As a first step, 
we define the Lie algebra elements 

Z ea ±e b = a ^ b em and Z 2ec = D+ e m, (3.65) 

where a,b,c E N ra and a < b. Observe that these matrices are labeled by the C n -type positive roots 
()A.2p . Now, we can formulate the main result of the paper. 

Theorem 4. The Lax matrix h2. 30\) of the classical hyperbolic BC n Sutherland model verifies the 
r-matrix Poisson bracket h3. 6°J\) with the q-dependent r-matrix 

r 12 (q) = 2 coth(a(g))X+> e ® - 2 £ . J Z a ® X~' 1 
' ^ sinn a o 

a (3.66) 
- ® £>+ + D c - £>") - ^(A+> £ ® X « ' e + X «' £ ® X « ' e )- 



2 We write Li = L ® 1, L 2 = 1 8> L, together with r 12 = £ r A - B T A ® T B , r 2 i = E r A - B T B ® T A , etc. 

3 For details on gauge transformations see e.g. [14], or Chapter 2 in [13] . 
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Proof. As dictated by (|3,64p . we inspect the formula of r\i (q, p) term-by-term. Using the basis given 
in the appendix, it is clear that the quadratic Casimir (|3.59p takes the form 

12 = ^2(D- ® D- - D+ ® D+) + Y.( X « e ® X « e ~ X a* ® X - ' e )- ( 3 - 67 ) 

Recalling (|2.27p and (|3.45|) . it is also clear that 

{^ l5 Y 2 } ext (n ) = Y, ((da R ) eQ (e Q T A ) + (dm) (eQ ^ E)) (e Q T A 0)) ® T A . (3.68) 



Now, from Lemma [1] it is immediate that 

^(da fi ) eQ (e Q T A ) ® T A = £ D+ ® D+ + £ X+< e (X+' e - coth(a(g))X-> e ) , (3.69) 

A c a,e 

whereas Lemma [2] leads to the formula 

Y,( dm \eQ^E))(e Q T A eO)®T A = -^^ (h c , sinh(ad )- 1 (T A ) n x) ® T A . (3.70) 

A Ac 

Remembering definition (|3.32j) . notice that the c-dependent part of H c can be rewritten as 

c— 1 n 

i(F c E* + EF*) = 2V2(D+ + X+f) + 2 ]T(X e + ;i ec + X e + ;} e J + 2 £ (X+'i^ + X+^J. (3.71) 

d=l d=c+l 

Therefore, upon introducing the Lie algebra element 

A 

we find easily that 

E(dm) (eQif(i;)) (e Q r A 0) ® T A = iljy ® + ]T — 1 Z a ® X"' 1 . (3.73) 

A a * 

By plugging formulae (|3.69|) and (|3.73p into (|3.68[) . the control over {<pi, l^ 2 } ext (uo) is complete. 
To proceed further, let us introduce the locally defined smooth functions 

a R (y, Y, p) = a R (y) and rh(y, Y, p) = m(y, p) ((y, Y, p) e P ext ). (3.74) 

Due to (|3.45p it is clear that (p = o R rh. Now, from (|2,27p we see that on ~p ext we have 

{(^)l,(^) 2 } ext =0, {(^) l5 (m) 2 } cxt ^O, {(m)i,(^) 2 } ext ^0, (3.75) 

therefore {(pi, <^ 2 } ext (-uo) = {(rh)i, (?Ti) 2 } ext (no) readily follows. Keeping our focus on this relationship, 
from (|2.25p and Lemma [2] it is immediate that 

V°Re(m c )(u ) = and V°Im(m c )(u ) = -S c /4. (3.76) 
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Note that the Lie algebra element S c (|3,32p can be represented as an appropriate commutator. Indeed, 
upon introducing the matrices 

1 - 

A c = - V(e c>(i - e d;C ) e u(n) and V c = diag(A c , A c ) £ t, (3.77) 

d=l 

we find immediately that H c = fi~ 1 [V c , £(E)]. Therefore, from ()2.27p and ()2.22p we obtain that 

{m c ,m d } cxt (u ) = {£(E), [V c ,V d ]) = 0, (3.78) 

for all c,d£ N n . Thus, we end up with the simple relationship {tpi, (f2} ext {uo) = 0. 

We conclude the proof with the observation that the term iljv ® ^(q) appearing in ()3.73p can be 
neglected, since it commutes with L\(q,p). Therefore, by simply plugging the above derived formulae 
into (|3.64p . the theorem follows. □ 

Switching to the standard basis {e^} of the matrix Lie algebra gl(N,C), the r-matrix (|3.66|) can 
be rewritten as 

n 

r i2{q) = ^ coth(g a - q b )(e a>b + e n+a , n+b ) (g> (e b>a - e n+b>n+a ) 

a,b=l 
n 



+ e n+a:b ) <8) (e 

n+b,a ^b,n+a) 



a,b=l 



(3.79) 



1 - 1 

+ + e b h + <S> (e a,b &n+a,n+b) 

2 smh(g a - q b ) 

a, o=L 

1 n 1 

— o / . . i / ; r(e a ,a + e n+a ,ri+a + e& b + e n+b , n+b ) (g) (e a n+b — e n+a b ) 

2 ^ i smh(g a + g 6 ) 

n 

+ ^ (e a ,6 <8> e n+b:n+a + e n+ajn+ 5 efc ia + e ajn+ fc ® eb >n+a + e n+a> b (8) e n+ t ia ) . 

a,6=l 

To conclude this subsection, notice that the above r-matrix can be seen as a generalization of the 
C n -type r-matrix constructed by Avan, Babelon and Talon. Indeed, up to a constant conjugation, the 
g-dependent part of (|3.79p can be identified with the r-matrix of the C n Sutherland model presented 
in [ID]. Nevertheless, as one can easily verify by inspecting the r-matrix Poisson bracket (|3.63[) . in the 
special case k = the ^-independent part of (|3.79|) can be safely neglected. In other words, with the 
specialization k = we can also recover the C n -type r-matrix of paper |10j . 



3.3 Lax representation of the dynamics 

Having constructed an r-matrix for the BC n Sutherland model, we can automatically provide a Lax 
representation for the dynamics as well. For this, we need the operator version of r\2(q) ()3.66p . which 
is defined via the natural identifications 

0®g = 0®0* = End(g). (3.80) 
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More precisely, the linear operator R(q) £ End(g) corresponding to the element 7*12(9) £ ® can De 
recovered from the formula 

R(q)Y = tr 2 (r 12 (q)Y 2 ) (Y G g), (3.81) 

where the linear operator tr2 defined by ta^X <g> Y) = ti(Y)X is the usual partial trace on the second 
factor. From (|3.66p and ()3.81|) it is immediate that 

R(q)Y = 2^coth(a(4))(X-' e ,Y)X+' e - 2 ]T — -^<X-\ Y~>Z Q - Y*. (3.82) 

Now, let us introduce the matrix- valued function B = i (L + -R£) defined on the phase space V s . 
Since the Lax matrix L (|2.30|) can be written as 



c=l c=l 

- 2/i ^ ( ^ e "- e 6 1 X e a +e b I 

l<^6<n \ Sinh ^ a ~ 96 ^ sinh (<?a + 9b) ) 

from (]3.82|) it follows that B has the block matrix structure 



(3.83) 



B 



S t 
T S 



(3.84) 



where S and T are appropriate u(n)-valued functions on V s . Namely, for their matrix entries we have 



. *^cosh(2g c ) + k _ . cosh(g a + q b ) _ . cosh(q a - q b ) 

1 . ; 2m \ ' a > b ~ 1 A i . ,0/ ', \ ' ^ a > b ~ J M . ,0/ 7' 

smh (2q c ) smh (q a + q b ) smh (q a - q b ) 



meanwhile 

. v + k cosh(2<7 c ) . /ll 



+ i/x V — -27 r + — r . ( 3 - 86 ) 



sinh (2g c ) ^ Vsmh(g c -g d ) smh (g c + 

where a,b,c G N n and a ^ b. Notice that i? is actually a t- valued map depending only on q. 

As we have discussed in Section [21 the reduced Hamiltonian corresponding to F 2 (|2.37p coincides 
with the Hamiltonian of the Sutherland model (|1.4[) . i.e. = (L, L) /4. By applying the Hamiltonian 
vector field X H s G X(V S ) on L, from the r-matrix Poisson bracket (|3.63p we obtain 

X H s[L} = ±[RL,L] = [B,L}. (3.87) 

That is, along each trajectory t \— > (q(t),p(t)) of the Sutherland dynamics the Lax equation 

L=[B,L] (3.88) 
holds. The above observation can be sharpened as follows. 

Theorem 5. A smooth curve in the phase space V s \1.2i) is an integral curve of the hyperbolic BC n 
Sutherland dynamics, if and only if, along the curve the Lax equation hS. 88\) is satisfied. 
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Proof. By applying repeatedly the identity 

cosh(x) 1 cosh(y) 1 1 / 1 1 



sinh 2 (x) sinh(y) sinh 2 (y) sinh(x) sinh(x + y) Vsinh 2 (x) sinh 2 (y)/ ' 
elementary algebraic manipulations lead to the formula 

cosh(g a + q b ) j 

I {Pa ~ Pb) . ,2, \^e a -e h T {Pa^Pb) 

l<a<b<n V Smh 

sinh 2 (2g c ) 2ec ^d 9c 



(3.89) 



rn rl o ft \ cos HQa ~ Qb) y -,i . , . . cosh(q a + q b ) 



(3.90) 



On the other hand, by differentiating L ()3.83p along an arbitrary smooth curve 

(?(*). P(*)) G ^ S with 

respect to time i, we find easily that 

4 - 2 ^ IS 5 £ „ («■ - + * + 

Hence, by comparing formulae (|3.90p and ()3.9ip . we conclude that the Lax equation is equivalent to 
the Hamiltonian equation of motion of the Sutherland model. □ 



4 Discussion 

Starting with the seminal paper [15] , a lot of effort has been devoted to explore the r-matrix structure of 
the Calogero-Moser-Sutherland many-particle systems. In this paper we contribute to this research 
area by constructing a dynamical r-matrix for the hyperbolic BC n Sutherland model with three 
independent parameters (jl .4|) . The outcome of our analysis is consistent with the results of [10] on 
the r-matrix structure of the C n Sutherland model with two independent coupling constants. We 
wish to mention that the authors of paper [16] have also constructed a dynamical r-matrix for a 
restricted class of the I?C n -type Sutherland models. More precisely, their results are valid under the 
same restriction on the coupling parameters that was sticked to these models in the fundamental work 
of Olshanetsky and Perelomov [I]. Recall also that the i?C n -type r-matrix in |16] explicitly depends 
on the coupling parameters. Note, however, that our r-matrix (|3.66p is independent of the coupling 
parameters, and so it is equally valid for the B n , C n and BC n Sutherland models, too. This 'universal' 
feature of ()3.66p naturally indicates a Yang-Baxter-type algebraic structure behind these models. We 
wish to investigate this important topic in future publications. A related open problem is to explore 
the relationship between our r-matrix and the i?C n -type Sutherland model with two types of particles 
(see e.g. [IT], M)- 

Similar questions arise in the context of the elliptic Calogero models, too. We have a fairly complete 
understanding of the r-matrix structure of Krichever's spectral parameter dependent Lax matrix [19 
for the A n -type model (see [20], [21] ). These elliptic r-matrices are dynamical objects, depending on 
the particle coordinates. As is known, one can even construct a Lax matrix for the elliptic A n -type 
model, which obeys an r-matrix Poisson bracket with Belavin's |22] non-dynamical elliptic r-matrix. 
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(The details on the elliptic case can be found in [23], whereas [M] contains an elementary account 
on the non-dynamical r-matrix structure of the degenerate A n -type models.) However, for the r- 
matrix structure of the elliptic BC n Calogero model only partial results are known. Namely, in [16J 
a dynamical r-matrix is constructed for the elliptic BC n Calogero model with the aforementioned 
restriction on coupling constants. For this restricted class of BC n models the Lax representation with 
non-dynamical r-matrix has been also investigated (see [25] ). Nevertheless, to our knowledge, the 
r-matrix of the most general three parameter dependent elliptic BC n Calogero model is still missing. 
It also appears to be an interesting open problem to provide r-matrices for the universal Lax operators 
constructed in [7J. An equally ambitious project would be to construct an r-matrix for Inozemtsev's 
[26J many-parameter dependent elliptic model, too. We hope that appropriate generalizations of our 
r-matrix (|3.66|) may play a role in clarifying these issues. 

To conclude the paper, let us recall that the Ruijsenaars-Schneider-van Diejen (RSvD) models 
(see e.g. [27], [28]) are natural generalizations of the Calogero-Moser-Sutherland (CMS) particle 
systems. The r-matrix structure of the ^4 n -type Ruijsenaars-Schneider models is well understood (for 
details on the elliptic models see e.g. [2JJ], [3D]), but for the generic non-^4 ra -type models even the Lax 
representation of the dynamics is missing. Quite surprisingly, the construction of a Lax matrix for 
the rational BC n RSvD model with three independent coupling parameters was carried out only in 
the recent paper [9]. Due to the dual reduction picture presented in [9], we expect that the r-matrix 
structure of the rational BC n RSvD model can be analyzed by the same techniques we outlined in 
Section [3j As for the ^4 n -type systems, it has been observed that in some sense the CMS and the 
RSvD models can be characterized by the same r-matrices (for details see [31]). It appears to be 
an interesting question whether the dual reduction picture behind the CMS and the RSvD models 
can provide a geometric explanation of this remarkable phenomenon. We wish to come back to these 
problems in future publications. 

A Convenient basis for u(n, n) 

As a supplementary material to the main text, in this appendix we present a convenient basis for the 
real Lie algebra q = u(n,n) adapted to the orthogonal decomposition (|2,10p . First, for each c 6 N n 
we define the linear functional 



can be seen as a family of positive roots of type C n . We also need the standard N x N elementary 
matrices e^/. Recall that for their matrix entries we have {ekijk' V = &k,k'°~i,i'- 
Now, for each c G N n we define the diagonal matrices 



e c :M n ^IR, q = (q u . . . , q n ) H- e c (q) = q c . 



(A.l) 



Clearly the set of functionals 



K+ = {e a ± e b | 1 < a < b < n} U {2e c | c G N n } 



(A.2) 



D+ = ^=(e c , c + 




c 




(A.3) 



Clearly {D^} is a basis in m, whereas {D c } is basis in a, 



satisfying the relations 




$c,d, (D c , D d ) = 5 Ctd . 



(A.4) 
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Next, for each c G N n we introduce the matrices 



-X-2e c — ^y^( e c,n+c ± ^n+c,c) ■ (A.5) 



Also, for all 1 < a < 6 < n we define the following matrices with purely real entries 

X e ^_ — ~F Cfc,a ^ &n+a,n+b ^n+b,n+a); 



^"e a -|_ ei) — 2 &b,n+a i ^-n+a,b ~F ^n+fc.a)] 



(A.6) 



together with the following matrices with purely imaginary entries 



-^e a — 2 (^ a >^ ^6, a &n+a,n+b ^n+b,n+a) i 



(A.7) 



Ag a + ei) 2 (&a,n+b "I - ^fe,n+a i Cn+a,6 i 6n+fe,a)- 

The set of vectors {Xq' € } forms a basis in m ± , meanwhile {Xa' e } is a basis in a^. Note that 

<A+' e ,X+/) = -<5 aia /<5 eie ,, <V- e , A"/) = 6 a ^6 e>e >. (A.8) 

Due to the orthogonality relations (|A.4|) and (|A.8|) . the construction of the corresponding dual basis 
is trivial. Keeping in mind the notation introduced in (|2.7p . it is worth mentioning that the above 
listed vectors satisfy the commutation relations 

[Q,Xt e ]=a(q)X^, (A.9) 

where q G R n , a G 7£ + and e G {r, i}. 
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